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Abstract 

We study the correlation functions in the branched polymer model. Although 
there are no correlations in the grand canonical ensemble, when looking at the 
canonical ensemble we find negative long range power like correlations. We pro- 
pose that a similar mechanism explains the shape of recently measured correlation 
functions in the elongated phase of 4d simplicial gravity[|^]. 

1 Introduction 

Recently some measurements of the curvature-curvature correlation functions in 4d sim- 
plicial gravity were presented [|^. The data seems to indicate negative and power like 
behavior both at the transition and in the elongated phase. The latter seems surprising 
as one would expect the exponential decay of correlations outside the critical region. Those 
issues are hard to study numerically and even harder analytically. To gain some insight 
into possible mechanism of those correlations we investigate here a simple geometrical 
model: branched polymers (BP). 

Besides its simplicity there are some other motivations for using this model. The BP 
describe the c — oo limit of 2d simplicial gravity H, . It is believed that BP also describe 
the elongated phase of 4d simplicial gravity [Q. However the exact correspondence be- 
tween the quantities measured in ref. and BP is not known. Nevertheless we can hope 
that the BP model captures the general features of the elongated phase of 4d simpUcial 
gravity. 



2 The model 

We consider here the ensemble of planar rooted planted trees. A tree is a graph without 
the closed loops. A rooted tree is a tree with one marked vertex (root) and planted tree 
is a tree whose root's degree (number of branches) is one. Two trees are considered as 
distinct if they cannot be mapped on each other by a continuous deformation of the plane 
(see fig. IH). We denote by n the total number of non-root vertices in the tree and by rij 
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Figure 1: Example of distinct planted rooted trees 
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Figure 2: Partition function 

the number of non-root vertices of the degree i. Then the (grand) partition function is 
defined as 

t)=Y: e-^y{T), p{T) = t^tT ■■■ (1) 
The T denotes the ensemble of all the trees. The equation for Z^n^ t) is (see fig. 0) 

oo 

z{ii,t) = e-^r{z{ii,t)) with r{z) = Y,uz'-^. (2) 

1=1 

The equation (pD can be rewritten as 

^(^) (3) 



and the critical point yUc corresponds to a value Zq where the right hand side of (|^) has a 
minimum. In the neighborhood of this point for a large class of parameters t the partition 
function behaves like 



Z(/i,t) ^ Zo(t) - Zi{t)^Afi + Z2{t)Ai2 + 0(A/i5) (4) 

with A/i = n — That is the only class of solutions of the equation (^ considered in 
this paper. Behaviour described by (§) is typical also for the elongated phase of the 4d 
simplicial gravity From the equation (Q) we can derive 
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Figure 3: Two-point function 

where ^"^^^ = 

3 The Two-Point functions 

First we consider a "volume-volume" correlation function [Q] 

G(/i,t;r)= Y: e-^-p{T) (6) 
T&rr 

where % is the ensemble of the trees with one point marked at distance r + 1 from the 
root. 

The smallest possible tree in %■ is a chain of r + 1 non-root vertices which we split 
into root, body and tail (see fig. ^). The weight of this chain is titgC"'-''''"^-*'^. All other 
configurations in % can be obtained from this chain by attaching trees in its non-root 
vertices(see fig. ^). Attaching n trees in the body or root part of the chain corresponds 
to a factor tn+2{n + ^)Z". The factor n + 1 counts the possible relative positions of the 
chain. Attaching n trees to the tail corresponds to a factor tn+iZ'"'. Finally our two-point 
function is 

^(/x, t; r) =e-'^(t2 + 2^2 + 3UZ^ + ■■■)■ 

■ e-('-i)''(t2 + 2t3Z + 3UZ^ + ■■■ y-^e~^{t^ + t^Z + ■ ■ ■ ) 

=e-^J^'{Z)e-^'''^^^J^'{Zf''-^^e-^'J^{Z) = (1 + Z 

Z' 

where J-''{Z) = and Z' = In the last equality we used the identity 

e~^^J-''{Z) = 1 + ^ obtained by differentiating the equation (|^) with respect to fi. 
We are interested in the Root-Tail correlation function defined as follows: 

i?T(/i, t;r) = J2 d{vi)d{vr+i)e-^y{T) (8) 

Tr 



(7) 
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where d{v) is the degree of vertex v. The form of this function can be easily obtained by 
modifying G{fi, t; r) 

RT{fi, t; r) =e-^(2 ■ ta + 3 ■ ^t^Z + 4 ■ 3^4^^ + ■■■)■ 

• e-('^-i)^(t2 + 2t3^ + 3t4^' + • ■ ■ y~^e~^{\ ■ h + 2 ■ + ■ ■ ■) 

1 g(^-(Z)Z^) ^, g(^(Z)Z) (9) 

Z dZ -A l^j e 

where we used the identity e~^F"{Z) = + -§n — again obtained by differentiating 
eq. (H) twice with respect to \i. 

Similarly we define two other functions 



Ay y2 7'^ 7" 7 

R{y.,t-T) = Y: d{v,)e-^y{T) = {2 + — + -,- — ^)(1 + -f-'^Z, 

TeTr Zj Zj Zj Zj 



(10) 



T(/., t; r) = 5: d{vr^ry-^y{T) = (1 + |^)'-(2Z + ^). (11) 



Let 



^ A(T)e-^>(r) 



(A)r = . (12) 

We define the normalized connected Root-Tail correlation function as 

RT,{^i,t;r) = {{d{vi) - {d{v^))rXd{vr+i) - (rf(t;.+i))rj)r. 

^ RT{^i,t;r) _ R{fi,t-r)T{fi,t;r) (13) 
G{fi,t;r) G{fi,t;ry 

It is easy to check that 

Rn{fx,t;r) = 0. (14) 

4 The canonical ensemble 

The functions defined in the preceding chapter are the discrete Laplace transforms of their 
canonical counterparts e.g. 

oo 

G'(/i,t;r) = ^e-'^''G(n,t;r) with G{n,t;r) = ^ p(T) (15) 

n=0 TeTr{n) 

where %{n) is the ensemble of the trees belonging to %. and having exactly n non-root 
vertices. To calculate the canonical functions from grand canonical we have to perform 
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the inverse of the discrete Laplace transform. This usually can not be done exactly and 
we proceed with a series of approximations. 

The first approximation is that of replacing the discrete transform (|1^) with the con- 
tinuous one. Then G{n, t; r) is given by the inverse Laplace transform which we calculate 
by the saddle point method. The details are presented in appendix A. Below we give 
results to the leading order in n (r <^ n). 

G{n, t; r) ^-^Zoxn'^ e''^" exp (-^x^] , (16) 
2^/iT \ An / 



2 ^ / Zl + Zf ZoiZl + 6ZoZ2-5Zf' 
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TT 



1 / Zo (Zo^ + 6Z0Z2 - 5Zf 
■ exp I \x 



(17) 



An \ Zl Z^ + Zi 



^ r. f Zf + Zl Zo(2Zl + 6Z0Z2 - AZl, , „ „ „ 
R(n, t; r) ^— fx ° - ^ ) n'^ e'^^" 

Zo {2Zl + 6Z0Z2 - AZl 
■^'^^^-Tn K'Yi 




T(n, t; r) ^1=Z. [x + |) n"! e-exp (-i- (x + |)^) (19) 



^i/ " V 4n V Zl. 

where we used x = ^^f ° 't^^ . 

The connected correlation function RTc(n,t;r) is defined as in (|T3|) 



RTciix,t;r) = {id{vi) - (rf(f i))r,{n))(c?(t^r+i) - (divr+i)) %{n) 

_ RT{n, t; r) _ R{n, t; r)T{n, t; r) (20) 
~ G{n,t;r) G(n,t;r)2 ' 

Before writing it down let us note that the terms proportional to (and so to r^) in the 
exponents cancel between numerators and denominators in RTc{n,t;r). For large n and 
r /n ^ 1 we can neglect the exponents. The resulting expression is 

RT in t- r) ^ ~ ^^'^ (21) 

' ' ^ zI{2tzI + zI) • ^ ' 

5 Examples and simulations 

Knowing the function J^{Z) it is easy to calculate (|21|). The equation (||) can be solved 
numerically if the analytic solution is not available. Below we give examples of two models. 
In both cases the correlations are negative. This was also the case for all other models 
we tested and we are persuaded that this is a general feature of the models with the 
expansion of the form (^) and all the weights t positive. 
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Figure 4: The results of MC simulations. 



1) ti — t, t2 — H — 1, ^4 — 



0, J^{Z) = t + Z + Z^, 
At 



RTc(n,t:r) = —- , _ 



(22) 



2) ti 



I 2 



r{Z) = |Lid(Z)§. For rf= 1, 



RTc{n; r) 



1.3029 



(0.725958r + 0.30219)2' 



(23) 



The formula ([21| ) is valid for n — >■ oo. To check what finite size effects are to be 
expected we performed the MC simulations of the second model {d = 1) with 256, 1024 
and 4096 non-root vertices. The results for the RTc are shown on fig. ^. We plotted 
the MC data, the large n predictions (formula (p3D) and the predictions for RTc without 
neglecting the exponents (the dotted lines). 



6 Discussion 

The appearance of correlations in the canonical ensemble is not surprising. What is more 
surprising is that those correlations do not vanish in the n oo limit. To understand 
better what is happening we calculate RTc{r = 1) for the first model from previous section 
with t = 1 by the explicit tree counting. In the fig. | we list all the relevant groups of trees 
together with their respective weights. The upper formulas refer to the grand canonical 
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ensemble and lower ones to canonical ensemble (valid for n^l). Note that the first two 
trees cannot appear in the canonical ensemble. If we calculate the RT^ using the grand 
canonical weights we get zero as expected. If we calculate the RT^ with grand canonical 
weights but exclude from the sum two first trees (those forbidden in canonical ensemble) 
we get a non-zero result depending on the value of w = exp(— /i). For the critical value 
w = I the result is — yI^-. Repeating the calculations with canonical weights we obtain 
— ^ which agrees with (^) for t = 1. This would indicate that the effect is due to the 
absence of small configurations in the canonical ensemble. 




Figure 5: Calculation of RTdr = 1) 



We have shown that a large class of BP models exhibits a long range negative power 
like correlations in the canonical ensemble. Those correlations are not finite size effects and 
survive in the n — > cxd limit. The shape of those correlations bears a striking resemblance 
to the shape of correlation functions measured in 4d simplicial gravity . Before making 
a detailed comparison one should keep in mind that the exact correspondence between 
the BP and 4d gravity is not known. The quantity measured in has only a qualitative 
resemblance to the quantity calculated here. We believe however that the same mechanism 
can be responsible for both. 



7 



7 Acknowledgments 



I would like to thank Zdzislaw Burda, Jerzy Jurkiewicz and Jan Smit for many helpful 
discussions and comments. This work is supported by the 'Stichting voor Fundamenteel 
Onderzoek der Matierie' (FOM). The numerical simulations were partially carried out on 
the PowerXplorer at SARA. 



A Calculation of G{n^ t] r] 

The inverse Laplace transform (continuous) of G{fi, t; r) is given by 

1 rc+ioD 

G(n,t;r)^— dftG{fi,t-r)e->^ (24) 

ZTTI J c—ioo 



where c > fic- The saddle point equation is 

d\og{G{fi,t;r)) 
(9/i 



-n. (25) 



Using the (^) and (|^) we can expand the left hand side of (|25|). Because we are interested 
in n ^ 1 we keep only the largest term of the expansion (in A/x) of the left hand side of 
(^Sf ). Solving this we obtain 



K^^^^z^^ (26) 

V? 



dHog{G{fis,t;r)) ^1 

y zi ^ 2Zoj 

Putting it all together we get the formula (|16|). Formulas for other two-point functions 
can be obtained in a similar way. 
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